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ABSTRACT 

We analyze PLANET and MACHO observations of MACHO 99-BLG-47, the first nearly-normal microlensing 
event for which high signal-to-noise-ratio data reveal a well-covered, short-duration anomaly. This anomaly occurs 
near the peak of the event. Short-duration anomalies near the peak of otherwise normal events are expected to arise 
both from extreme-separation (either very close or very wide), roughly equal-mass binary lenses, and from planetary 
systems. We show that the lens of MACHO 99-BLG-47 is in fact an extreme-separation binary, not a planetary 
system, thus demonstrating for the first time that these two important classes of events can be distinguished in 
practice. However, we find that the wide-binary and close-binary lens solutions fit the data equally well, and cannot 
be distinguished even at Ax^ = 1. This degeneracy is qualitatively much more severe than the one identified for 
MACHO 98-SMC-l because the present degeneracy spans two rather than one dimension in the magnification field 
and does not require significantly different blending fractions. In the appendix, we explore this result, and show that 
it is related to the symmetry in the lens equation. 

Subject headings: binaries: general - gravitational lensing - planetary system 



1. Introduction 

The hallmark of planetary microlensing events is a short deviation from an otherwise normal, point-source/point-lens (here- 
after PSPL) event. Mao & Paczyriski (1991) showed that extrasolar planets could be detected from such events, and Gould & 
Loeb (1992) gave an explicit prescription for how the planet/star mass ratio q (^ 1) and the angular separation d (in units of 
the angular Einstein radius 6'e) could be reconstructed by decomposing the event light curve into its "normal" and "perturbed" 
components. 

Work during the ensuing decade has elucidated many additional subtleties of planetary light curves, but their fundamental 
characterization as briefly perturbed PSPL events has remained intact. Of particular importance in the present context, Griest & 
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Safizadeh (1998) showed that events with small impact parameter (uq ^ 1; where uq is the minimum separation between the 
source and the lens center of mass in units of 9e) probe the so-called "central caustic" of the lens geometry, making them much 
more sensitive to the presence of planets than the larger impact-parameter events analyzed by Gould & Loeb (1992), which probe 
the outer "planetary caustic". 

These central-caustic events are of exceptional importance, even though they are intrinsically rare. They are rare simply 
because the central-caustic is much smaller than the planetary caustic, so the great majority of planet-induced deviations (of 
fixed fractional amplitude) are due to planetary caustics. However, the probability of detecting a planet is much greater in small 
impact-parameter events, partly because the source is guaranteed to pass close to the central caustic, which almost by definition 
is near the center of the lens geometry (u = 0) and partly because even the sensitivity of planetary caustics is enhanced for 
Mo <C 1. (Here u denotes the source position on the sky, normalized to 6'e, with respect to the lens center of mass.) By contrast, 
higher impact-parameter events miss the central caustic, and they are likely to miss the planetary caustic as well because it 
lies in a random position relative to the source trajectory. Because of their higher sensitivity to planets, and because they can 
be recognized in real time, low impact-parameter events are monitored more intensively than typical events by microlensing 
follow-up networks, which in turn further enhances their sensitivity to planets. 

Central-caustic events, like their planetary-caustic cousins, involve a short deviation from an otherwise normal PSPL light 
curve. The major difference between these two classes of planetary events is that central-caustic anomalies always occur near 
the peak, whereas planetary-caustic perturbations can occur anywhere on the light curve, and are typically expected on the wings 
of the light curve. Of particular importance, for central-caustic events, there is no simple prescription for extracting d and q by 
decomposing the light curve into "normal" and "perturbed" components and it is unclear to what degree these two parameters are 
degenerate. 

Another important, albeit accidental, discovery was that planets could give rise to perturbations that are not short compared 
to the event timescale. In the course of their search for planetary perturbations among 43 approximately PSPL events, Gaudi et 
al. (2002, see also Albrow et al. 2001b) found one event, OGLE-1999-BUL-36, that was asymmetric in a way that was consistent 
with the presence of a planet. They argued, however, that the asymmetry was also consistent with parallax effects induced by 
the Earth's motion around the Sun of the type analyzed by Gould, Miralda-Escude, & BahcaU (1994), and that in general it 
would be extremely difficult to distinguish between the two possible causes of such an asymmetry. They concluded that, in 
most cases, microlensing searches are not able to distinguish between parallax and a weak, asymmetric planetary perturbation, 
and consequently, all such "detections" should be ignored. This reduces the sensitivity of microlensing searches to planets, but 
only by an extremely small amount since, as Gaudi et al. (2002) showed, long-timescale asyrmnetric perturbations account for 
less than ~ 1% of all planetary events. Hence, the long-timescale asymmetric events also confirm in a way the basic paradigm: 
planetary perturbations have short durations relative to the parent light curve, and in the rare cases for which they do not, they are 
not recognizable as a planetary perturbation anyway. 

However, not all short timescale deviations are due to planets and therefore the mere detection of such an anomaly does not 
prove the presence of a planet. Gaudi & Gould (1997) showed that close binaries (d <^ 1) give rise to light curves that are virtually 
identical to PSPL events, except when the source comes very close to the lens center of mass (|m| <C 1). Hence, for events with 
Uq ^ 1, the light curve looks "normal" except for a brief deviation near the peak. Qualitatively, this is exactly the same behavior 
as that of central-caustic planetary events. Similarly, light curves of wide-binary (d » 1) events can also take the same form if 
one - and only one - of the caustics lies very close to the source's passage. Indeed, a close correspondence between a certain 
pair of close-binary and wide-binary events was discovered both theoretically (Dominik 1999) and observationally (Albrow et al. 
1999; Afonso et al. 2000). It remains an open question under what conditions these various types of events can be distinguished 
from one another. If central-caustic planetary events could not be distinguished from close- and/or wide-binary events, it would 
seriously undermine planet searches in high-magnification events and hence would call into question the basic strategy adopted 
by microlensing follow-up groups (e.g., Albrow et al. 1998). 

Here we analyze the light curve of the microlensing event MACHO 99-BLG-47, the first intensively monitored event with a 
short-lived, high signal-to-noise-ratio deviation from an otherwise normal PSPL light curve. We identify the lens as an extreme- 
separation binary rather than a planetary system, thereby showing that, at least in this case, the two classes can be clearly 
distinguished. We also show that both wide- and close-binary solutions fit the data equally well, implying that although planetary 
perturbations can be distinguished from those arising from extreme- separation binaries, the discrimination between very close 
and very wide binaries may be difficult in practice. 
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2. MACHO 99-BLG-47 

The initial alert for the microlensing event MACHO 99-BLG-47 was issued by the MACHO collaboration on 1999 July 
23.'^ The event was located about 18° from the Galactic center along the disk (I = 17°59'10", b = — 1°57'36") and reported 
to be rather faint at basehne (V = 21.5, R = 20.3). PLANET began monitoring MACHO 99-BLG-47 right after the electronic 
alert with the expectation that it would be a very high magnification event. PLANET detected anomalous behavior in the Ught 
curve during the first week of August, and subsequently issued an anomaly alert on 1999 August 4.^^ 



2.1. Observations and Data 

The PLANET light curve of MACHO 99-BLG-47 consists of observations made from three different southern sites: the 
Elizabeth 1 m at the South African Astronomical Observatory (SAAO), Sutherland, South Africa; the Canopus 1 m near Hobart, 
Tasmania, Austraha; and the Yale/AURA/Lisbon/OSU (YALO) 1 m at Cerro Tololo Inter-American Observatory (CTIO), La 
Serena, Chile. The observations were carried out in two bands: / (at all three sites) and V (at SAAO only). The event was most 
intensely monitored in 1999 August during and following the photometric peak and the anomaly, but we obtained data through 
mid-September and also when the event was close to baseline early in the 2000 season as the source came out from behind the 
Sun. After the usual reductions, we perform photometry by two independent methods: a direct fit to the point spread function 
(PSF) using DoPHOT (Schechter, Mateo, & Saha 1993) and difference imaging using ISIS (Alard & Lupton 1998; Alard 2000). 
For details, see also Albrow et al. (2000a). Because of the different characteristics of the two methods and varying frame qualities, 
we recover a different number of photometric measurements. Before any elaborate efforts to clean and select the data, we have 
303 points from DoPHOT photometry and 297 from ISIS photometry (Table 1) for the 1999 observations, as well as 8 points 
observed from YALO during early 2000 season (DoPHOT photometry only). 

The photometric observations taken at SAAO are calibrated to the standard Johnson-Cousins system with respect to E-region 
stars (Menzies et al. 1989) observed contemporaneously with MACHO 99-BLG-47 at SAAO on 1999 July 31. The instrumental 
magnitudes of observations from Canopus and YALO are not calibrated expUcitiy, but rather we allow independent source fluxes 
for different sites and bands in our modehng, which provide the relative photometric scalings among them. 

We also include in our analysis the pubUcly available MACHO photometry, which comprises 77 Vmacho and 90 i?MACHO 
points from the 1999 season as well as 213 Vmacho and 227 i?MACHO points taken between 1995 and 1998. The latter were 
used to constrain the baseline and to check for source variability. We find no evidence for such variability. 

The difference photometry is placed on an absolute scale by a linear regression from the PSF-based flux (DoPHOT) to 
the differential flux (ISIS), after removing obvious faulty points from the two data sets based on the reports by each reduc- 
tion/photometry package. We find that, except for the Canopus data, the residuals of the PSF-based flux from the regression line 
are strongly correlated with the size of the seeing disk, but that most of this correlation is removed by adding a linear seeing 
correction. For example, for the SAAO /-band observations, we detect the presence of a non-zero linear seeing dependence term 
of ~ l.TFbase arcsec"^ (or ^ lOFg arcsec"^), with the signal to noise ratio being as large as 30. Here Fg and .Fbase are the 
net flux of the lensed source star alone and the total flux in the same PSF when there is no magnification. 

In order to minimize the effects of known systematics, we optimize the data set by rejecting outUers and reevaluating the 
sizes of the photometric error bars. We also include a correction term for the correlation between the seeing and the blended 
flux that enters the same PSF with the lensed source when we fit the observed flux to a specific magnification model. These 
procedures are described in detail and fuUy justified in several earlier PLANET papers (Albrow et al. 2000a,b, 2001a,b; Gaudi et 
al. 2002). Briefly, we first construct an initial clean subset of the data by rejecting 3-(j outliers with respect to the linear regression 
from DoPHOT to ISIS flux, and then determine a reference model based on this subset. Once we have a reference model, we 
include all the data and follow the same iterative procedure of outlier removal and error reseating described in Albrow et al. 
(2001a). The final data set used for the subsequent analysis reported in the present paper contains 276 DoPHOT-reduced and 266 
ISIS-reduced PLANET points as well as 166 MACHO points (one Vmacho point rejected), all from the 1999 season (Table 1). 
The data obtained during seasons other than 1999 (flie 2000 season YALO /, and flie 1995 - 1998 season MACHO Vmacho and 



'^ftp://darkstar.astro.washington.edu/macho/Alert/99-BLG-47/ 
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^MACHO bands) are included in the analysis as a combined single point for each set. In doing so, 4 data points (all -Rmacho in 
the 1998 season) are eliminated as outliers. 



3. Analysis and Model 

The light curve of MACHO 99-BLG-47 (Fig. 1) is that of a normal high-magnification PSPL event with a short-lived (^ 

3 days) deviation near the peak, of the type predicted for planets (Griest & Safizadeh 1998) or binaries with extreme separations. 
We therefore develop a method for probing the space of lens geometries of this type and then search for the minimum within 
this space. 

We begin by excising the anomalous points near the peak and fitting the remaining light curve to the degenerate form of the 
Paczynski (1986) profile that obtains in the high-magnification limit (Gould 1996), 



F{t) = -—^^^=— + i^base • (1) 

Here F{t) is the instantaneous flux, Fpcak is the peak flux above the baseline of the model, Fbase is the baseline flux, to is the time 
of the peak, and teff is the effective width of the light curve. These degenerate parameters are related to the standard parameters 
by -Fpeak = Fg/uo, -Fbase = + Fh, and ^eff = uqIe, where is the Einstein timescale, i.e., the time required for the source to 
move an Einstein radius. We find to ~ 2451393.6 (Heliocentric Julian Date), tes ~ 1 day, and that -Fpeak roughly corresponds 
to /peak ^ 15-3 (Fig. [1]). However, the exact values of those parameters are quite dependent on which points we excise. Note 
that though the baseline magnitude is measured to be /base = 19.1, there is at this point essentially no information of how the 
baseline flux divides into source flux Fg and blended flux F^, which is why it is necessary to fit the light curve to the degenerate 
profile (eq. [1]). 

Next, we reinsert the excised points near the peak and note that the maximum (over-magnified) deviation from the degenerate 
fit is ^ 0.5 mag (Fig. 1). We then establish a set of initial trial event geometries as follows. For each diamond-shaped caustic 
produced by various geometries of {d,q) pairs, we examine the magnification as a function of distance from the "caustic center" 
along each of the three directions defined by the four cusps of the caustic (one direction is redundant due to the reflection 
synmietry with respect to the binary axis). Here for computational simplicity, we use an analytic proxy point for the "caustic 
center," defined as follows. For close binaries, we adopt the center of mass of the lens system for the caustic center On the 
other hand, for wide binaries, the position of one component of the binary shifted toward the other component by [d{l + q~^)]~^ 
is chosen as the caustic center. [These choices can be understood from the approximation developed in Appendix A. See also 
Dominik (1999).] Then, at each point ti,, the point of the cusp-axis crossing, we take the ratio of the actual magnification to 
that predicted for the corresponding PSPL lens: a point mass lens located at the caustic center with the same mass as either the 
whole binary (close binaries) or the nearest lens alone (wide binaries). We proceed with the search only if this ratio lies in the 
interval [1.2,1.8]. We then choose a source trajectory perpendicular to the cusp axis as an initial trial model, with uq = \u, \ and 
tE = tes/\u,\, and use four-parameter (a, ^e, to, and uo) downhill simplex (Press et al. 1992) to search for a best-fit model to 
minimize x^. which is determined by a linear- flux fit to the model of lens magnifications (and thus the corresponding Fg and F^ 
are found immediately), subject to the constraint of fixed id,q). Figure 2 shows the result of this search (based on ISIS solutions) 
in a contour plot of x^-surface over the (d,g')-space. The x^ is rising around all the boundaries shown in Figure 2, except toward 
lower q. While the x^-surface flattens with Ax^ — 32 as the mass ratio q becomes very small (q < 0.01) for d ~ 0.065 and 
d ~ 15, we find no evidence of a decrease of the x^ as q becomes smaller than 0.01. Rather it asymptotically approaches 
Ax^ ^ 32. 

We find well-localized minima of x^ over the searched id,q)-spacc, one of close binaries (d ~ 0.13) and the other of wide 
binaries (d ~ 11.3), whose exact parameters depend slightly on whether we use the ISIS or DoPHOT photometry (Tables 2 and 
3). We adopt the ISIS solutions in the subsequent discussion (and they are also what is shown in Fig. 2) because the light curve 
shows less scatter and consequently, the errors for the model parameter determinations are smaller. Despite the combination of 
the dense coverage near the peak by the PLANET data and the extensive baseline coverage by the MACHO data, the final two 
ISIS models (i.e., the close binary and the wide binary) are essentially indistinguishable: Ax^ = 0.6 for 412 degrees of freedom, 
with the close binary having the lower x^ of the two. In Appendix A, we discuss this degeneracy in further detail. 

Finally, we also note that uotE and F^/uo from the best-fit binary-lens model parameters are not quite the same as the tcff and 
-Ppeak derived from the initial degenerate form of the light curve fit. These discrepancies are traceable to small differences between 



-5- 



a true Paczynski (1986) curve and its degenerate approximation (eq. [1]), and are not due to differences between the binary and 
corresponding PSPL event, which are very small except around the peak. Since the parameters derived from equation (1) function 
only as seeds for simplex, and since the final x^-surface is very well-behaved, these discrepancies in input values do not influence 
the final result. 



4. Discussion 

The mass ratios of the best fit models are g = 0.340 ± 0.041 (close binary) and q = 0.751 ± 0.193 (wide binary), which 
are well away from the regime of planetary companions {q < 0.01). By comparison, the ratio of the duration of the anomalous 
portion of the light curve {t 

anom ^ days) to the Einstein timescale (^e 160 days)^^ is r — ^anom/^E 0.02. For light curves 
perturbed by planetary caustics, one typically finds q ^ (Gould & Loeb 1992). This relation clearly does not apply to the 

caustics of extreme-separation binaries. 

Models with q < 0.01 (and q > 100) are formally rejected at Ax^ — 32, which is significant but not in itself an overwhelm- 
ing rejection of the planetary hypothesis. Hence, we examine the least model for q = 0.01 (Ax^ = 32.3) for its plausibility 
and the origin of statistical discriminating power. We discover that the difference of x^ is mostly from the MACHO data between 
HJD 2451310 and 2451360, approximately 2 months prior to the photometric peak. However, we also find that the "planetary" 
model exhibits an extremely high peak magnification (Amax ~ 15000), and consequently, requires the event to be much longer 
(fg — 40 yr) and more extremely blended (4 ~ 25.7, Fg/Fbase — 0.2%) than the already unusually long and highly blended 
best-fit models. In addition, there is a clear trend of increasing peak magnification (and thus timescale and blend as weU) as q is 
lowered beyond 0.01. This follows from the fact that the observed timescale of the anomaly in the light curve essentially fixes 
the source movement relative to the caustic. However, as q becomes smaller, the size of the caustic relative to the Einstein ring 
shrinks, and therefore, the timescale of the event, that is the time required for the source to cross the Einstein ring, increases. 
(Note that this behavior causes a mild continuous degeneracy between q and the blending.) These parameters determined for 
the "planetary" model are extremely contrived and highly improbable a priori. Furthermore, the timescale associated with the 
"planet" component, tp = q^^'^t^ ~ 4 yr is much longer than that of typical stellar lenses, which further reduces the plausibility 
of the planetary interpretation. In summary, while in simple statistical terms, the star/planet scenario is not overwhelmingly 
disfavored relative to extreme- separation binaries, we nevertheless can conclude that it is highly unlikely that this event is due to 
a star/planet system. 

From the standpoint of refining microlensing planet detection strategies, it is important to ask how one could have discrim- 
inated between the planetary and extreme- separation binary solutions with greater statistical significance. As noted above, most 
of the discriminating power came from the MACHO data points on the rising wing of the light curve, even though (or in a sense, 
because) these had the largest errors and the lowest density of the non-baseline coverage. That is, the precision PLANET pho- 
tometry over the peak and falling wing "predicts" the rising wing for each of the models, but the noisier MACHO data can only 
roughly discriminate between these predictions. Hence, the key would have been to get better data on the rising part of the light 
curve. In practice this is difficult: the MACHO data are noisier exactly because they are survey data, and one does not know to 
monitor an event intensively until the light curve has actuaUy started to rise. 

The long duration of the event, — 160 days (close binary) or t-E — 220 days (wide binary), may lead one to expect 
that the event would show some sign of parallax effects (Gould 1992; Mao 1999; Soszyfiski, et al. 2001; Bennett et al. 2001). 
Similarly, the close passage of the source to a cusp could in principle give rise to finite source effects (Gould 1994; Nemiroff & 
Wickramasinghe 1994; Witt & Mao 1994), as for example was recently observed in the case of EROS BLG-2000-5 (An et al. 

2002). Indeed, the combination of parallax and finite-source effects permitted An et al. (2002) to measure the mass of a microlens 
for the first time. We have therefore searched for both parallax and finite source effects, but find no significant detection of either. 

Finally, we examine the position of the source on the color-magnitude diagram (CMD; Fig. 3). The most prominent feature 
found in the CMD is a track of stars running diagonaUy in the same direction as would the main sequence. Since the field is in 
the Galactic disk, this feature is probably not a true main sequence, but rather likely to be a "reddening sequence", that is, an 
ensemble of mostly disk tumoff stars at progressively greater distances and correspondingly greater reddenings. The baseline 



While the value here is for the close binary, the value for the wide binary is essentially the same because the relevant mass of the corresponding point-mass 
lens for the wide binary is not the combined mass but the mass of the single component that the source passes by. 
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"star" (combined light of the source and blended light) lies within this sequence toward its faint/red end although its position is 
seeing-dependent because of the seeing-correction term. The position indicated in the CMD is plotted assuming median seeing. 
However, the source itself [/g = 20.9, {V — I)s = 1.94] lies 2 mag below the baseline "star," in a region of the CMD that 
is well below the threshold of detection. If the baseline were mainly composed of source light, and the "blended light" were 
simply source light that had been falsely attributed to blending by a wrong model, then the real source would lie within the well 
populated "reddening sequence", and the timescale would be much shorter, ~ 30 days. We therefore searched for solutions 
with little or no blending. However, we find that these are ruled out with Ax^ = 2440. Furthermore, the best-fit models with no 
blended light are still extreme-separation binaries and not planetary systems. Since we have only a crude understanding of the 
CMD, no strong argument can be made that the position of the source is implausible. If, in fact, the observed track of stars in 
Figure 3 is really the reddening sequence of disk tumoff stars, the source position is consistent with a low-mass main-sequence 
star lying > 2 mag below the tumoff. This would also explain the lack of finite source effect. Moreover, we note that in disk 
fields, long events are not at all uncommon (Derue et al. 2001) because the observer, source, and lens are all moving with roughly 
the same velocity. 

Regarding the large amount of blended light required to fit the observed light curve, we note that the highly significant 
seeing effect in the PSF-based photometry (see § 2. 1) is also evidence that the event is strongly blended. In addition to the direct 
confirmation of the strong seeing effect from the comparison between the differential flux and the PSF-based flux measurements, 
we independently detect a significant seeing correction when we fit the observed flux to the magnification model. For SAAO 
/-band observation, the seeing correction determined from the model fit using DoPHOT flux is slightly smaller (~ QFg arcsec"^) 
than the value derived from the regression between DoPHOT and ISIS flux. The seeing correction derived by fitting ISIS data to 
the model is basicaUy consistent with zero. 



5. Binary Lens vs. Binary Source 

Multiple-peak events like the one seen in Figure 1 can be caused by a binary source (Griest & Hu 1992) rather than a 
binary lens. In general, one expects that such events will be chromatic, since the colors of the two sources will not usually 
be the same. By contrast, MACHO 99-BLG-47 is achromatic: the difference in color of the two wings of the light curve is 
AV — I = 0.01 ± 0.05. In the limit that one component of the binary is completely dark, the event will be achromatic, and 
yet can have still have multiple peaks caused by rotation of the binary during the event. In this case, however, there will be a 
series of roughly equally spaced peaks that gradually die out as the event declines (Han & Gould 1997), contrary to the distinct 
double-peaked behavior seen in Figure 1. Hence, a binary-source explanation appears a priori implausible. 

Nevertheless, it is possible in principle that the source has two components of nearly identical color. We search for binary- 
source models, with either static or slowly moving components, but find that the best fit model has Ax^ = 1167, which clearly 
rules out such models. The basic problem is that the observed second peak is very sharp given its height and the timescale of 
the declining light curve. Binary-source light curves that fit these latter features tend to have a width that is closer to that of the 
dotted curve in Figure 1. Recall that this curve represents a single-lens degenerate fit to the non-peak data. We conclude that the 
explanation for the double-peaked behavior of the light curve is that it is a binary-lens rather than a binary-source event. 



6. Conclusion 

MACHO 99-BLG-47 is the first microlensing event with a short-lived, high signal-to-noise-ratio anomaly, characteristics 
that could betray the existence of a planet around the lensing star. Nevertheless, we conclude that the lens of MACHO 99-BLG-47 
is not a planetary system, but an extreme-separation (very close or very wide) binary composed of components of similar mass, 
based on the result of the light curve fit as well as the extreme value of event duration and blending fraction required for any 
plausible "planetary" fit. 
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A. Close and Wide Binary Correspondence 

Dominik (1999) noticed that the Chang & Refsdal (1979, 1984) lens (CRL; see also Schneider, Ehlers, & Falco 1992) well 
approximates the binary lens system in the vicinity of each individual lens component of extreme wide- separation (d 1) 
binaries and in the vicinity of the secondary of extreme close-separation (d <C 1) binaries while the quadrupole lens (QL) ap- 
proximation works nicely near the center of mass of extreme close-separation binaries. He further showed that the size and shape 
of the caustics and the behavior of the magnification field associated with them are very similar among the true binary lenses, the 
CRL approximation, and the QL approximation if the shear, 7, of the CRL and the absolute eigenvalue of the quadrupole moment 
tensor, Q (for simplicity, hereafter Q will be just referred to "the quadrupole moment"), of the QL are both very small and their 
numerical values are close to each other (7 ~ Q ^ 1). These findings also imply the existence of a correspondence between 
a wide binary with 7 — d^^qu,{l + Qw)^^ and a close binary with Q = dfqc{l + Qc)^'^, sometimes referred to as a d ^ d^^ 
correspondence. (Here and throughout this appendix, we use the subscript c and ^ to distinguish the parameters associated with 
the close and wide binaries.) Here we rederive the basic result for this correspondence, and explore it more thoroughly. 

The general form of the binary lens mapping equation can be expressed in notation utilizing complex numbers (Witt 1990), 

C = ^ - - - — — , (Al) 

Z — Zi Z — Z2 

and the magnification associated with a single image can be found by, 

-1 



A= 1 



(A2) 



where C,, z, zi, and Z2 are the positions of the source, the image, and the two lens components normalized by the Einstein radius 
of some mass, and ei and 62 are the masses of lens components in units of the same mass. For the extreme close-binary case, 
if one chooses the center of mass as the origin and the binary axis as the real axis and sets the combined mass to be unity, then 
zi = dc€2, Z2 = —dc€i, ei = (1 -|- qc)~^, and €2 = 9c(l + 9c)~^, and the lens equation (Al) becomes 

C = Z-- — 1 — (A3) 

z — dce2 z + dcEi 

« + + (A4) 

Z Z'^ z* 

Here the monopole term (z~^) is the same as for a PSPL and the first non-PSPL term is the quadrupole (z^'^), which acts as the 
main perturbation to the PSPL if the quadrupole moment is small; Q = d^eie2 — d^qdl + qc)^"^ ^ 1- For a given source 
position Q, let the image position for the PSPL be zq. Then, under the perturbative approach, the image for the QL approximation 
is found at ^; = 2:0 -h 5zc, where \5zc\ <C |-2o| ^ 1. Using the fact that zq is the image position corresponding to C, for the PSPL 
(i.e., C, = zq — Zq^), one obtains 
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from equation (A4) and taking only terms linear in 5zc and Q. Here, we note that Q ^ 0(df ) so that the second term in the 
bracket is lower order than Q^. In order to find the magnification for this image, we differentiate equation (A4), and find. 



dC 

dz 



1 

~2 



3Q 



1 



4(1 - qc) 4 
3(1 + qc) z 



Then substituting z = zq + 5zc, for which 5zc is given by equation (A5), into equation (A6), we obtain. 



dz 



+ Q 



3ko|4-2|2o| 



1 / 1 



^0 



^0 



4|zo|*-2|zo| 



^0 



1 \ 1 



l + qc 



(A6) 



(A7) 



Note that the total magnification for the given source position is usually dominated by one or two images found close to the 
critical curve. Thus, we consider only the case for which the non-perturbed PSPL images lie close to the unit circle, so we have 
1^0 1 = 1 + A and A <C 1. Then, we find the expression for the inverse magnification for the QL approximation (up to the order 

ofdl). 



4A-2Q(1 + 1 

-^0 



30(4 + 4 

^0 ^0 



1 



i^oi - 1) - O3?(^o"') + Irrfl'^'^^^^^""'^ 



(A8) 



For the extreme wide-binary case, one can rewrite the lens equation (Al) as 



. 1 

z z + di 

^ 1 qw ^ Qw 

z di dl 



q-w 2 , 
dl^ 



{dw > 1^1) • 



(A9) 
(AlO) 



Here the position and the mass of the first lens component are the origin and the unit mass so that zi = 0, Z2 = —di, ei = 1, 
£2 = qw, and di = {1 + qwY^^di^. We note that, apart from the constant translation, the first non-PSPL term here is essentially 
the shear, 7 — qwdi^ = qwd~'^{l + qw)~^ for the CRL approximation. Analogous to the extreme close binary, if 7 <C 1, the 
image position for the CRL approximation can be found by the perturbative approach, but here the corresponding PSPL source 
position would be C + q-wd^^ = zq — Zq^. Then, the image deviation Sz^ = z — ZQof CRL from PSPL is 



5zu 



7 1 



^0 



4) {l+q^Y'^du 



Using this result and the derivative of equation (AlO), 



1 



dz^ z''^^ 



we find 

dz 



+ 7 



1^0 I 



■2|zop-3 



^0 



1 



(i + g^)V2 4 



2|0o|6 + 2|zo|4-4|zo| 

ko|'*-i 



(l + g,,)i/2rf,^ 



(All) 



(A12) 



(A13) 



From the same argument used for the QL approximation of the extreme close binary, we can set |2:o| = 1 -h A, and then the 
inverse magnification for the CRL approximation (up to the order of d^"*) is 



1 



1 



1 



4A-27 3 + ^ +37 -3+^ 



(|zo|-l)-73?(^o"') + 



1 



2{l + q^Y/^du 



[l + qn,Y/Hu 



zo') 



(A14) 



By comparing equations (AS) and (A14), we therefore establish the magnification correspondence (up to the order of d^ ~ d^^) 
between the close binary with Q = d!^qc{l + qc)~'^ and the wide binary with 7 = d~^g^(l -|- qw)~^ when Q — 7 ^ 1. 

For the two PLANET ISIS models for MACHO 99-BLG-47, we obtain Q = 3.40 x 10"^ for the close binary solution 
{dc, qc) = (0.134, 0.340) and 7 = 3.35 x 10~^ for the wide binary solution {dw,qw) = (11-31, 0.751). Hence, the observed 
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degeneracy between the two PLANET models (see § 3 and Tables 2 and 3) is the clearest example of this type of correspondence 
observed so far. While Afonso et al. (2000) reported that two different binary lens models, one of a close binary and the other of a 
wide binary, can fit the observed light curve of the (caustic-crossing) binary lens event, MACHO 98-SMC-l, one can infer from 
their figure 8, that the degeneracy between the two models exists only for the specific light curves (which is essentially a particular 
one dimensional slice of the magnification field over the source plane) but not for the magnification field in the neighborhood 
of the caustic as a whole. This is obvious from the relative rotation of the two caustics. In addition, the source magnitudes for 
the two models of MACHO 98-SMC-l differ by ~ 0.18 mag. By contrast, the difference of the predicted 4 between the two 
models of MACHO 99-BLG-47 is only - 0.02 mag. In fact, we get Q = 6.5 x 10"^ (the close binary) and 7 = 1.8 x 10"^ 
(the wide binary) for the two models of MACHO 98-SMC-l, and thus, although the degeneracy of the MACHO 98-SMC-l Ught 
curve is somehow related to the correspondence of d <-> d~^,it cannot be completely explained simply by the argument in this 
appendix, and it should be investigated further for its origin. On the other hand, the degeneracy of MACHO 99-BLG-47 is the 
first definitive observational case of the correspondence between extreme separation binaries. This can be also seen in Figure 4, 
which illustrates the similarity between the magnification fields for the two models of MACHO 99-BLG-47. 

The very low Ax^ between the two solutions despite the excellent data implies that it is extraordinarily difficult to break 
this degeneracy with photometric data. Gould & Han (2000) showed that for MACHO 98-SMC-l the two solutions were also 
astrometrically degenerate, at least for data streams lying within a few of the peak. They argued that this astrometric degen- 
eracy, like the corresponding photometric degeneracy, was rooted in the lens equation. However, as shown in this appendix, the 
correspondence between the equations describing close- and wide-binaries is purely local (See also Dominik 2001). For example, 
there is a constant-offset term in equation (AlO), which does not give rise to any local photometric or astrometric effects, but 
which must "disappear" at late times. Hence there must be a late-time astrometric shift between the two solutions. Such a shift 
was noted explicitly by Gould & Han (2000) for the case of MACHO 98-SMC-l, and these are likely to be a generic feature of 
close/wide corresponding pairs of solutions. 

We also plot the lines of id,q) pairs that have the same shear 7 as the best-fit wide-binary model or the same quadrupole 
moment Q as the best-fit close-binary model on x^-surface contour plot shown in Figure 2. While the iso-shear line for wide 
binaries lies nearly parallel to the direction of the principal conjugate near the best-fit model, it is clear from the figure that the 
condition of Q — 7 -C 1 alone does not define the observed well-defined two-fold degeneracy, which involves the additional 
correspondence between higher order terms beyond the quadrupole moment (~ d^) and the pure shear d^^). Further compar- 
ison between equations (A8) and (A14) indicates that there exists a magnification correspondence up to the order of ~ d~^ 
if the condition dc{l — (Zc)(l + Ic)^^ = d^,^{l + qw)^^^"^ is also satisfied in addition to Q = 7 ^ 1. We find that, for the two 
PLANET models, 4(1 - gc)(l + Sc)"^ = 6.6 x 10"^ (the close binary) and ^-^(1 -h g,^)"^/^ = 6.7 x IQ-^ (the wide binary). 
Hence, we conclude that, in fact, these two conditions, 

dldl{l + q^) = ^{l+qcf, (A15) 
ded^il + q^y/^ = , (A16) 

1 - gc 

define a unique correspondence between two extreme separation binaries. We also note that the images (not shown) of the iso- 
Ax^ contours for the close binary models of MACHO 99-BLG-47 under the mapping defined by the above two relations follows 
extremely closely the corresponding iso-Ax^ contour for wide binary models, except for the difference of Ax^ — 0.6 offset 
between two solutions. 



'^Albrow et al. (2001a) also found a similar degeneracy while they modeled the Ught curve of OGLE-1999-BUL-23, but they could discriminate the two 
models to better than Ax^ ~ 128. 
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Table 1 . Photometric data of MACHO 99-BLG-47 



telescope filter # of total points # of points analyzed median seeing error rescaling factor 

(arcsec) 



PLANET (DoPHOT) 



Elizabeth 1 m 


I 


98 


88 


1.81 


1.19 




V 


26 


24 


1.79 


1.02 


Canopus 1 m 


I 


51 


42 


2.71 


1.43 


YALO Im 


I 


128 


122 


1.96 


0.772 


PLANET (ISIS) 


Ehzabeth 1 m 


I 


105 


98 


1.65 


1.35 




V 


21 


20 


1.96 


1.83 


Canopus 1 m 


I 


57 


38 


2.59 


1.75 


YALO 1 m 


I 


114 


110 


1.83 


0.919 


MACHO 


Mt.Stromlo 50" 


-Rmacho 


90 


90 




1.04 




Vmacho 


11 


76 




0.709 'yO.730'' 



"^with respect to PLANET/DoPHOT model 
''with respect to PLANET/ISIS model 
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Table 2. PLANET close-binary model of MACHO 99-BLG-47 





PLANET flSIS") + MACHO 


PLANET fDoPHOTl + MACHO 


d 


0.134 ±0.009 


0.121 ±0.009 


Q 


0.340 ± 0.041 


0.374 ± 0.058 




163 ± 26 days 


183 ± 32 days 


to 


1393.9331 ±0.0071 


1393.9309 ±0.0075 




(8.6 ± 1.1) X 10'^ 


(7.5 ± 1.3) X 10-3 




294?99 ± 0?25 


294?45 ± 0?25 




412.13 


420.45 


dof 


412 


423 



Note. — We note that the fact that ~ dof results from our rescaling of the photometric errorbars. However, we use the 
same scaling factor here and for models described in Table 3 so the indistinguishability between two models is not affected by 
this rescaling. The uncertainties for parameters are "1-cr error bars" determined by a quadratic fit of surface. 

'^the closest approach to the binary center of mass 
^the HeUocentric JuUan Date - 2450000. 

'^the binary center of mass lying on the right hand side of the moving source 



Table 3 . PLANET wide-binary model of MACHO 99-BLG-47 



parameters 


PLANET (ISIS) + MACHO 


PLANET (DoPHOT) + MACHO 


d 


11.31 ±0.96 


12.79 ± 1.05 


q 


0.751 ±0.193 


0.917 ±0.288 


tE" 


220 ± 37 days 


253 ± 47 days 


to 


1393.9113 ±0.0072 


1393.9133 ±0.0074 




(6.5 ±1.0) X 10-3 


(5.50 ±0.92) X 10-3 


ad 


295?46 ± 0?24 


294?83 ± 0?24 




412.74 


420.50 


dof 


412 


423 



"^with respect to the Einstein radius of the combined mass 

''The closest approach to the caustic center. See § 3 for the definition of the caustic center, 
'^the HeUocentric JuUan Date - 2450000. 

^the caustic center lying on the right hand side of the moving source 
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Fig. 1. — Observations and models of MACHO 99-BLG-47, all scaled to calibrated Cousins /-band with blending as registered 
by SAAO DoPHOT photometry. Shown are PLANET /-band from SAAO, Canopus, YALO, PLANET V-hzad from SAAO, 
and MACHO Vmacho and /?macho- The ISIS reduced PLANET points are plotted after applying the transformation to the 
absolute scale derived from the linear regression between difference photometry and PSF-based photometry. The solid curve 
shows the final close-binary lens model fit (Table 2) to the data while the dotted curve shows the "best" degenerate form of 
point-source/point-lens light curve (eq. [1]) fit to a high-magnification subset of the data that excludes the anomalous points near 
the peak. The half-magnitude offset between this curve and the data is the main observational input into the algorithm to search 
for the final model (see § 3). Note that on the scale of the figure, the wide-binary solution is completely indistinguishable from 
the close-binary solution, i.e. the solid curve can represent both the close-binary and the wide-binary lens models equally well. 
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0.05 0.1 0.5 1 5 10 

d 

Fig. 2. — Contour plot of x^-surface over (d,q) space based on solutions for PLANET (ISIS) and MACHO data. The binary 
separation d is in units of the Einstein radius of the combined mass, and the mass ratio q is the ratio of the farther component to 
the closer component to the source trajectory (i.e., g > 1 means that the source passes by the secondary). Contours shown are of 
Ax^ = 1, 4, 9, 16, 25, 36 (with respect to the global minimum). We find two well-isolated minima of x^, one in the close-binary 
region, ((i,(7)=(0. 134,0.340), and the other in the wide binary-region, ((i,(7)=(l 1.31,0.751) with Ax^ = 0.6 and the close binary 
being the lower x^ solution. Also drawn are the curves of models with the same Q as the best-fit close-binary model and the 
same 7 as the best- fit wide-binary model (see Appendix A). 
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Fig. 3.— CMD of 3' X 3' field surrounding MACHO 99-BLG-47 (I = 17°99, b = -1?96) derived from SAAO observations. 
Shown are the positions of the baseline "star" (S+B) and the lensed source (S) in the absence of lensing (close-binary model). 
The position of the lensed source in the wide-binary model differs from this by substantially less than the size of the symbol. 
The errors are also smaller than the symbols. The majority of stars in the CMD are most likely turn-off stars seen at increasingly 
greater distances in the Galactic plane, and hence at correspondingly greater reddenings. 
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Fig. 4. — Grayscale plot of the difference of the normalized flux fields, 2\Fi^ — Fc\/{Fs^u,A^ + Fg^c^c) of the two PLANET 
models around the caustic. The two fields are scaled and aligned so that the source trajectories coincide with each other The 
resulting shape and the relative size of the caustics are remarkably close to each other [c.f., fig. 8 of Afonso et al. (2000) and fig. 6 
of Albrow et al. (2001a)]. Also drawn are the contours of zero difference (dotted line) and 5% difference (solid line). While the 
actual source trajectory (y = 0) naturally traces well the zero-difference line, the flux difference would also be extremely smaU 
for other trajectories through the region, except very near the caustic. 



This figure "fig4.jpg" is available in "jpg" format from: 
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